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Multifractal detrended cross-correlation analysis for two nonstationary signals 
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It is ubiquitous in natural and social sciences that two variables, recorded temporally or spatially in 
a complex system, are cross-correlated and possess multifractal features. We propose a new method 
called multifractal detrended cross-correlation analysis (MF-DXA) to investigate the multifractal 
behaviors in the power-law cross-correlations between two records in one or higher dimensions. The 
method is validated with cross-correlated ID and 2D binomial measures and multifractal random 
walks. Application to two financial time series is also illustrated. 
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Fractals and multifractals are ubiquitous in natural 
and social sciences [l[. The most usual records of ob- 
servable quantities in real world are in the form of time 
series and their fractal and multifractal properties have 
been extensively investigated. There are many methods 
proposed for this purpose 0, Q • For a single nonstation- 
ary time series, the detrended fluctuation analysis (DFA) 
can be adopted to explore its long-range autocorrelations 

0, [H and multifractal features [y] . The DFA method can 
also be extended to investigate higher-dimensional fractal 
and multifractal measures Q. 

There are many situations that several variables are 
simultaneously recorded that exhibit long-range depen- 
dence or multifractal nature, such as the velocity, temper- 
ature and concentration fields in turbulent flows [a, [9 , 1C | , 
topographic indices and crop yield in agronomy [111 . 1 12 1 , 
asset prices, indexes, and trading volumes in financial 
markets [13, [l4|. Recently, a first method named de- 
trended cross-correlation analysis (DXA) is proposed to 
investigate the long-range cross-correlations between two 
nonstationary time series, which is a generalization of the 
DFA method [ill. Here we show that the DXA method 
can be generalized to unveil the multifractal features 
of two cross-correlated signals and higher-dimensional 
multifractal measures. The validity and potential util- 
ity of the multifractal detrended cross-correlation anal- 
ysis (MF-DXA) method is illustrated using one- and 
two-dimensional binomial measures, multifractal random 
walks (MRWs), and financial prices. 

Consider two time series {xi} and {yi}, where i = 

1, 2, ■ • • , M. Without loss of generality, we can assume 
that these two time series have zero means. Each time 
series is covered with M s = [M/s] non-overlapping boxes 
of size s. The profile within the vth box [l v + 1, l v + s], 
where l v = (v — l)s, are determined to be X v (k) = 
£j=i x{l v +j) and Y v {k) = J2j=i y( l v+j), k = 1, • • • , s. 
Assume that the local trends of {X v (k)} and {Y v (k)} are 
{X v (k)} and {Y v (k)}, respectively. There are many dif- 
ferent methods for the determination of X v and Y v . The 



trend functions could be polynomials 0. The detrend- 
ing procedure can also be carried out nonpar ametri cally 
based on the empirical mode decomposition method [lq . 
The detrended covariancc of each box is calculated as 
follows 
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Y v (k)-Y v (k) (1) 



The qth order detrended covariance is calculated as fol- 
lows 



F xy (q,s)= -$>«( S ) 
v—1 

when q and 

F xy (0,s) = exp 
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We then expect the following scaling relation 

F xy (q,s)~s h *«M . (4) 

When X = Y, the above method reduces to the classic 
multifractal DFA. 

In order to test the validity of the MF-DXA method, 
we construct two binomial measures from the p-modcl 
with known analytic multifractal properties as a first ex- 
ample jig . Each multifractal signal is obtained in an 
iterative way. We start with the zeroth iteration g = 0, 
where the data set z(i) consists of one value, = 
1. In the gth iteration, the data set {z^(i) : i = 
1,2,- ■■ ,23} is obtained from z^(2fc + l) = pz^- 1 ) [2k + 
l)and^(2fc) = {l-p)z^- 1) (2k) for k = 1,2,- ■■ ,2 s " 1 . 
When g — > oo, z^ 9 '(i) approaches to a binomial measure, 
whose scaling exponent function h zz (q) has an analytic 
form [H [H 



H zz (q) = l/q-\og 2 [p q + (l-p) q ]/q 



(5) 
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In our simulation, we have performed g = 17 iterations 
with p x = 0.3 for x(i) and p y — 0.4 for y(i). The cross- 
correlation coefficient is 0.82. We find that F xy , F xx 
and F yy all scale with respect to s as power laws. Note 
that there are evident log-periodic oscillations decorat- 
ing power laws, which is an inherent feature of the con- 
structed binomial measures The best estimates of 
the power-law exponents are obtained when s is sampled 
log-periodically [2fjj . The resultant power-law exponents 
h xy , h xx and h yy are illustrated in Fig.[TJ The MF-DFA 
analysis gives h x 
also find that 



10' 



, : {q) = H xx {q) and h yy {q) = H yy {q). We 



hx V {q) = [h xx {q) + h yy (q)}/2 



(6) 



For monofractal ARFIMA signals, this relation with q 
2 is also observed Il5| . 
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FIG. 2: (color online). Multifractal nature in the power-law 
cross-correlations of two MRWs. (a) Power-law scaling in F xy , 
F xx and F yy with respect to s for q = 2 and 5; (b) Power-law 
exponents h xy , h xx and h yy . 



FIG. 1: Scaling exponents h xy (q) estimated using the mul- 
tifractal detrended cross-correlation analysis of two cross- 
correlated binomial measures generated from the p-model. 
The numerically estimated exponents h xx (q) and h yy (q) ob- 
tained from the multifractal detrended fluctuation analysis of 
x(i) and y(i) are also presented, which well match the analyt- 
ical curves H xx (q) and H yy (q). This example illustrates the 
relation h xy (q) = [h xx (q) + h yy (q)} /2. 

As a second example, we consider the multifractal ran- 
dom walks (MRWs) [I]} . The increments of an MRW are 
e[k]e^ k \ where e[k] and u>[k] are uncorrelated and u>[k] is 
a white noise. In order to generated two cross-correlated 
MRWs, we can generate two time series e x and e y pos- 
sessing the properties in the MRW formwork and rear- 
range e y such that the rearranged series e y has the same 
rank ordering as e x (22|. We generate two MRW signals 
of size 2 16 with A 2 = 0.02 for x{i) and A 2 = 0.04 for 
y(i), whose cross-correlation coefficient is 0.70. When q 
is negative, no evident power-law scaling is observed for 
Fxy(s), which has great fluctuations. When q is positive, 
nice power- law scaling is observed for F xy , F xx and F yy , 



as illustrated in Fig.[2^a) for q = 2 and 5. The power-law 
exponents h xy , h xx and h yy are illustrated in Fig. [2jb). 
We see that Eq. (|6|) holds in repeated numerical exper- 
iments. However, this relation does not hold for some 
other realizations of MRWs. 



We now apply the MF-DXA method to the daily clos- 
ing prices of DJIA and NASDAQ indexes. For compar- 
ison, we have used the same data sets and same scaling 
range as in Ref. 15[. No evident power-law scaling is 
observed for negative q values. For positive q, we see 
power-law dependence of F xy , F xx and F yy against time 
lag s. Two examples are illustrated in Fig. EJa) for q = 2 
and 5, where the case of q = 2 reproduces the results in 



Ref. [1 51 ] - The power-law exponents h xy , h xx and h yy are 
depicted in Fig. [2](b) , which are nonlinear functions with 
respect to q. We see that each time series of the absolute 
returns possesses multifractal nature and their power-law 
cross-correlations also exhibit multifractal nature. 

We can generalize the ID MF-DFA to the 2D version 
and its extension to higher dimensions is straightforward. 
Consider two self-similar (or self-affine) surfaces of iden- 
tical sizes, which can be denoted by two arrays x(i,j) 
and y(i,j), where i = 1,2, • • • , M, and j = 1, 2, • • • ,7V. 
The surfaces are partitioned into M s x N s disjoint square 
segments of the same size s x s, where M s = [M/s] 
and N s = [N/s]. Each segment can be denoted by 
Xv,w or y v<w such that x VtW (i,j) = x(l v + i,l w + j) and 

Vv,w{i,j) = y(lv +i,lw + j) for 1 < i,3 < s , where 
l v = (v — l)s and l w = (w — l)s. 

For each segment x v>w identified by v and w, the cu- 
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FIG. 3: (color online). Multifractal nature in the power-law 
cross-correlations of the absolute values of daily price changes 
for DJJIA and NASDAQ indexes in the period from July 1993 
to November 2003. (a) Power- law scaling in F xy , F xx and 
F yy with respect to s for q — 2 and 5. The scaling range 
is the same as in Ref. [l^]. (b) Dependence of the power- 
law exponents h xy , h xx and h yy as nonlinear functions of q, 
indicating the presence of multifractality. There is no clear 
relation between these exponents. 



mulative sum X ViW (i,j) is calculated as follows: 

i j 

Xv,w(i,j) = X! X! x v,w(ki,ki) , (7) 
fc 1= i k 2 =i 

where 1 ^ i, j ^ s. The cumulative sum Y VtW (i,j) can be 
calculated similarly from y VjW . The dctrended covariancc 
of the two segments can be determined as follows, 

- s s 

F v , w (s) = ^3 EE [ Xv > w -^v,w{i,j) 

■ (8) 



1=1 j=l 

X 



where X VjW and are the local trends of X v>w and 
Y v>Wl respectively. The trend function is pre-chosen in 
different function forms The simplest function could 
be a plane u(i,j) = ai + bj + c, which is adopted to 
test the validation of the method. The overall dctrended 
cross-correlation is calculated by averaging over all the 



segments, that is, 

Fxy(q,s) 



^ s s v =lw=l ) 



1/9 



(9) 



where q can take any real value except for q = 0. When 
q = 0, we have 



F xy (q, s) = exp 



«^*-"» ■ < 10 ' 



according to L'Hospital's rule. The scaling relation be- 
tween the detrended fluctuation function F xy (q,s) and 
the size scale s can be determined as follows 



F xy {q, s) 



,{q) 



(11) 



Since N and M need not be a multiple of the segment 
size s, two orthogonal strips at the end of the profile 
may remain. Taking these ending parts of the surface 
into consideration, the same partitioning procedure can 
be repeated starting from the other three corners [23| ■ 

It is noteworthy to point out that the order of cumula- 
tive summation and partitioning is crucial in the analy- 
sis of two- or higher-dimensional multifractals. Consider 
the point located at (l v + i, l w + j) in the box identified 
by v and w, where 1 $5 i,j ^ s. The cumulative sum 
X(l v + i,l w + j) can be expressed as follows 

lv lw 

+ j) = X 

E E x ( k i,k 2 ) 

kl=l fe 2 =l 

l v l w +j l v +i l w 

+ E E ^i> fc 2)+ E £*(*!,*%,). (12) 



For any pair of (z, j), X v>w (i,j) is localized to the seg- 
ment a;„ iUJ , while X(l v + i, l w + j) contains extra informa- 
tion outside the segment as shown above, which is not 
constant for different i and j and thus can not be re- 
moved by the detrending procedure. We find that the 
power-law scaling is absent if X(l v + i,l w + j) is used 
in Eq. (J7J). This observation is analogous to the case of 
higher-dimensional detrended fluctuation analysis Q ■ 

We now present numerical experiments validating the 
two-dimensional multifractal dctrended cross-correlation 
analysis. There exist several methods for the synthesis 
of two-dimensional multifractal measures or multifractal 



rough surfaces [2J|. The most classic method follows a 
multiplicative cascading process, which can be either de- 
terministic or stochastic QJ, 25, 2(| 27 1. The simplest one 
is the p model proposed to mimic the kinetic energy dis- 
sipation field in fully developed turbulence [17[ . Starting 
from a square, one partitions it into four sub-squares of 
the same size and assigns four given proportions of mea- 
sure pn, pi2, P21, and P22 to them. Then each sub-square 
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is divided into four smaller squares and the measure is re- 
distributed in the same way. This procedure is repeated 
g times and we generate multifractal "surfaces" of size 
2 s x 2 9 . The analytical expression of H xx (q) or H yy (q) 
for individual multifractals is 

H(q) = [2 - log 2 (pfi + P \ 2 + P q 21 +Ph)} /q ■ (13) 

In our simulation, we have usedpn = 0.10, p± 2 = 0.20, 
p 2 i = 0.30, and p 22 = 0.40 for X and p n = 0.05, 
p 12 = 0.15, p 21 = 0.20, and p 22 = 0.60 for Y. We 
find that the cross-correlation coefficient between the two 
multifractals depends linearly on the generation number 
g: c = —0.0408,g + 0.9528, where the value of R-squared 
is 0.9997. The 95% confidence intervals for the slope and 
intercept arc [-0.0415, -0.0402] and [0.9489,0.9566], re- 
spectively. In our numerical experiment, we have used 
g = 12, which gives c = 0.48. Very nice power-law be- 
haviors are confirmed in F xy (q, s), F xx (q, s), and F yy (q, s) 
with respect to s for different values of q. The resul- 
tant power-law exponents h xy (q), h xx (q) and h yy (q) are 
illustrated in Fig. 0] marked with open circles, squares 
and triangles, respectively. We find that the relation 
h xy (q) = [h xx {q) + h yy (q)] /2 holds. 
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FIG. 4: Multifractal detrended cross-correlation analysis of 
two cross-correlated synthetic binomial measures from the p- 
model. The size of each multifractal is 4096 x 4096 and the 
cross-correlation coefficient is 0.48. The numerical exponents 
h xx (q) and h yy (q) obtained from the multifractal detrended 
fluctuation analysis of X and Y locate approximately on the 
analytical curves H xx (q) and H yy (q). This example illustrates 
the relation h xy (q) = [h xx (q) + h yy (q)\ /2. 

In summary, we have proposed a multifractal de- 
trended cross-correlation analysis to explore the mul- 
tifractal behaviors in power-law cross-correlations be- 
tween two simultaneously recorded time series or higher- 
dimensional signals. The MF-DXA method is a com- 
bination of multifractal analysis and detrended cross- 
correlation analysis. Potential fields of application in- 
clude turbulence, financial markets, ecology, physiology, 
geophysics, and so on. 
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